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We discuss the relation between the Polyakov loop and the chiral order parameter at finite tem- 
perature by using the Gocksch-Ogilvie model with fundamental or adjoint quarks. The model is 
based on the double expansion of strong coupling and large dimensionality on the lattice. In an 
analytic way with the mean field approximation employed, we show that the confined phase must be 
accompanied by the spontaneous breaking of the chiral symmetry for both fundamental and adjoint 
quarks. Then we proceed to numerical analysis to look into the coupled dynamics of the Polyakov 
loop and the chiral order parameter. In the case of fundamental quarks, the pseudo-critical tem- 
perature inferred from the Polyakov loop behavior turns out to coincide with the pseudo-critical 
temperature of the chiral phase transition. We discuss the physical implication of the coincidence 
of the pseudo-critical temperatures in two extreme cases; one is the deconfinement dominance and 
the other is the chiral dominance. As for adjoint quarks, the deconfinement transition of first order 
, persists and the chiral phase transition occurs distinctly at higher temperature than the deconfine- 

ment transition does. The present model study gives us a plausible picture to understand the results 
from the lattice QCD and aQCD simulations. 
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"st 1 I. INTRODUCTION 
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Recent developments in the heavy ion collision experiment are going to enable us to reach phase transitions that 
■ the early universe must have undergone, i.e., the phase transition from hadronic matter to a quark-gluon plasma 
(T) ' (QGP) 1]. From the theoretical point of view what is known as the QGP phase transition can consist of two distinct 
phase transitions. One is the color deconfinement transition and the other is the chiral symmetry restoration or the 
chiral phase transition. The color liberation associated with deconfinement is a phenomenon distinguishable from 
the chiral phase transition in principle, but it is likely that the mechanism of confinement would be closely related 
to the chiral dynamics. According to several theoretical studies 0, in fact, confinement should be accompanied by 
Oh! the spontaneous chiral symmetry breaking. In other words, the critical temperature relevant to the deconfinement 
transition (denoted by Td) is inevitably lower than that relevant to the chiral phase transition (denoted by T x ), that 
' is; Td < T x , if two phase transitions are to be distinguished. 

Actually it is difficult to shed light on underlying physics in the relation between two phase transitions. It is partly 
because most of the effective models involving the chiral dynamics have shortcomings in describing confinement 
physics, and partly because the criterion to characterize the deconfinement transition cannot go together with the 
realization of the chiral symmetry. We can define an order parameter for the deconfinement transition only when 



X. 

there is no dynamical quark in a theory, or quarks are infinitely massive (m q — > oo). In contrast to it, the chiral 
symmetry is manifested when quarks can be regarded as massless {m q ~ 0). Thus two phase transitions must lie in 
the opposite limits in regard to the current quark mass. 

In the absence of dynamical quarks pure gluonic systems have the center symmetry and the Polyakov loop, L, 
provides a criterion for confinement 0, The expectation value of the Polyakov loop traced in the fundamental 
representation can be related to the free energy f q with a static fundamental quark placed in a thermal medium, 
that is; (Tr c L) = exp[—/3f q ], where Tr c denotes the trace over the color space and (3 is the inverse temperature [38|. 
Under the center transformation the traced Polyakov loop, Tr c L, is twisted by an element of the center group (i.e. 
the Z(N C ) group for an SU(iV c ) gauge theory). In the confined phase where f q — oo, the center symmetry is realized 
to result in (Tr c T) = 0, while the spontaneous breaking of the center symmetry occurs in the deconfincd phase where 
the expectation value of the Polyakov loop takes a finite value corresponding to f q < oo. 

For the purpose of looking into the spontaneous breaking of the center symmetry, the effective action in terms of the 
Polyakov loop has been pursued by the strong coupling expansion as well as by the perturbative calculation. Indeed, 
it has been revealed that the center symmetry is spontaneously broken in a perturbative regime at high temperature, 
which is actually embodied in what is often called the Weiss potential 0. As the temperature is lowered, unstable 
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modes can develop to make the expectation value of the Polyakov loop diminished |8|. Unfortunately, however, we 
cannot reach an appropriate description of the deconfinement transition in any perturbative way because the phase 
transition is a non-perturbative phenomenon and moreover the perturbative calculation would have no information on 
the Haar measure 8, 9J. The Haar measure in the context of confinement physics means a sort of the Faddeev-Popov 
determinant jlCj . or precisely speaking, the Jacobian of the variable transformation from the gauge potential to the 
Polyakov loop. As stressed also in Ref. the nature of the group integration contained in the Haar measure can 
be responsible for confinement and, in fact, it is the case in the strong coupling expansion on the lattice. It has been 
found that the resultant effective action in the lowest order of the strong coupling expansion leads to a second order 
phase transition for the SU(2) gauge theory ^1 an d a nrs t order phase transition for the SU(iV c > 3) gauge theories 
|l3j. These results in the strong coupling expansion are consistent with the anticipation derived from the analogy to 
spin systems based on the center symmetry [3( . They are in agreement with the results from the lattice simulation as 
well [4J. 

When massless quarks in the fundamental representation are present in a theory, the center symmetry is broken 
explicitly and the Polyakov loop would no longer serve as an order parameter to identify the deconfinement transition 
|14j . Then, another symmetry, i.e., the chiral symmetry plays an important role in hadronic properties. There are 
many effective approaches based on the chiral symmetry such as the linear-sigma model, the Nambu-Jona-Lasinio 
model, the chiral random matrix model and so on |l5j . The spontaneous breaking of the chiral symmetry with the 
staggered fermion on the lattice has already been investigated by means of the double expansion of strong coupling 
and large dimensionality [16l Il7| . The results are surprisingly acceptable not only qualitatively but quantitatively 
also. 

For a moderate value of the quark mass, neither the deconfinement transition nor the chiral phase transition has 
any well-defined order parameter. In the lattice QCD simulation, it has been observed that the Polyakov loop and the 
chiral order parameter show crossover behaviors in the presence of dynamical quarks 0, as if the magnetization in 
spin systems under an external magnetic field. To be interesting, the lattice QCD results suggest that the simultaneous 
transitions for deconfinement and chiral restoration take place at the same pseudo-critical temperature, that is; the 
respective peaks of the susceptibility for the Polyakov loop and for the chiral order parameter are found at the same 
temperature. 

In Ref. [I9|, Satz emphasized that such coincidence of the pseudo-critical temperatures can be understood in the 
following way: At low temperature where the chiral symmetry is spontaneously broken, the explicit breaking of the 
center symmetry is suppressed by the heavier mass of constituent quarks rather than the lighter mass of current 
quarks. Consequently the expectation value of the Polyakov loop stays small in the confined or chiral broken phase 
at low temperature. Once the constituent quark mass drops off in the chiral symmetric phase at high temperature, 
the explicit breaking of the center symmetry is enlarged to result in a finite value of the Polyakov loop. Although this 
qualitative argument is developed further in Ref. |20l |. any quantitative demonstration had not been given until our 
analysis on the Gocksch-Ogilvie model pll |22|] . It should be noticed that, according to Satz's argument, the behavior 
of the Polyakov loop would be sensitive to the chiral phase transition rather than the deconfinement transition. 
In order to disclose the underlying relation between deconfinement and chiral restoration, some legitimate order 
parameters are indispensable to characterize the deconfinement transition even in the presence of dynamical quarks. 
However, defining an order parameter for the deconfinement transition is a quite difficult question because of the 
thermal excitation of dynamical quarks, particularly under the thermodynamic limit at finite temperature 23]. 

Alternatively, it is interesting to consider the system with dynamical quarks in the adjoint representation since 
adjoint quarks would preserve the center symmetry. Actually the color structure of adjoint quarks is identical to 
that of gluons. Thus the Polyakov loop serves as an appropriate order parameter for the deconfinement transition 
with respect to fundamental quarks in a similar way as in the pure gluonic theories. The chiral phase transition for 
dynamical adjoint quarks and the deconfinement transition for static fundamental quarks can be defined separately in 
this case. In the lattice simulation with adjoint quarks (lattice aQCD), it has been found that the critical temperature 
for the chiral restoration is distinctly higher than that for the deconfinement transition [24[ . In this paper we apply 
the Gocksch-Ogilvie model to investigate the system with adjoint quarks as well. 

This paper is organized as follows: We explain what the Gocksch-Ogilvie model is with emphasis on its physical 
meaning in Sec. [H] We give computational details in Add. 1X1 additionally. In Sec. IIIII we elaborate the mean field 
approximation employed in the present analysis. Then in an analytic way we investigate the effects on the chiral 
dynamics arising from the Polyakov loop behaviors. Sec. IIVI is devoted to presenting the numerical results and 
discussing their physical implications. Finally, we conclude our model study in Sec. 
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FIG. 1: Intuitive pictures for the effective action of the Gocksch-Ogilvie model. The torus consists of the thermal S 1 and the 
spatial R . The arrowed line represents the link variable. The filled and blank circles are the quark and anti-quark excitations, 
(a) In the spatial directions a quark and an anti-quark propagate as a meson as a result of the group integration with respect 
to the spatial link variables, (b) Thermal excitations of dynamical quarks are allowed by the non-vanishing Polyakov loop at 
finite temperature. 



II. GOCKSCH-OGILVIE MODEL 

The Gocksch-Ogilvie model is given by the effective action 0,H2; (see App. 1X1 for detail derivations) 
S$[L,\] = -J^Tr c L(n)Tr c Lt(^) + £ X(n)V(n, m)A(m) 

n.n. m,n 

~ T ^ ^ ln [ C0Sh ( N r E ) + \ ( L(r) + LWt )] . (1) 
ft 

where L(n) is the Polyakov loop in the fundamental representation defined on the lattice by 

L{n) = ]J U d {ti,n d ) (2) 

n d =a 

in the d dimensional space-time. denotes the Polyakov loop in the r representation (r = fund, for the fundamental 
representation and r — adj. for the adjoint one). The Polyakov loop in the adjoint representation can be expressed in 
terms of the fundamental ones as 



4f ° = 2tr 



(3) 



where the matrices t a form a fundamental representation of the SU(iV c ) group and the indices a, b run from 1 to 
iV c 2 — 1. The normalization is determined to satisfy tr[i a ib] = 5 a b/2. 

In the effective action X(n) is the meson field, Tr c the trace with respect to the color indices, and Nf the number 
of flavors which is taken as Nf = 2 throughout this paper, dim(r) is the number of dimension in the r representation 
of the SU(iV c ) group, that is; dim(fund.) = N c and dim(adj.) = N% ~ 1. The temperature T, the meson hopping 
propagator V(n, m), the strength of the nearest neighbor interaction J, and the quasi-quark energy (constituent quark 
mass) E are defined respectively as 



i 

N~a ' v " , "" / ~ 2{d 



n,m-\-j n,m—j/ 
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J = e~^ T , E = si ^-UJLJ: x + mq ) (4 ) 



with j running from 1 to d — 1 only in the spatial directions, a stands for the string tension fixed at the empirical 
value at zero temperature, i.e., a = (425 Me V) 2 . The lattice spacing a and the current quark mass m q are treated as 
the model parameters. 

Physical meaning for each term in the effective action Q is straightforward to understand. The first term is the 
gluonic contribution in terms of the Polyakov loop. An order-disorder phase transition described by this gluonic action 
signifies the transition from the center symmetric confined phase to the symmetry broken deconfined phase. The next 
term involving A(n) is the meson interaction term in the spatial direction resulting from the group integration with 
respect to the spatial link variables as shown in Fig. ^ (a). The last logarithmic term corresponds to the thermal 
excitation of dynamical quarks as depicted intuitively in Fig. ^ (b). This interpretation will become tangible if once 
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the Polyakov loop is set to be unity by hand. The logarithmic term then takes a familiar form of the free energy 
for fcrmionic ideal gases. 

In the case of fundamental quarks the center symmetry is explicitly broken by the thermal quark excitation with 
fractional baryon number. Actually the expansion of the logarithmic quark contribution in the power of the Polyakov 
loop is essentially equivalent to the hopping parameter expansion pEj and then the power of the Polyakov loop counts 
the winding number of quark paths along the thermal torus, namely, the quark number in the system. The leading 
contribution in such an expansion in terms of the Polyakov loop is ~ Tr c Z/ cosh(iV T .E), which can be regarded as 
the symmetry breaking term under an external magnetic-like field in analogy to spin systems. The magnitude of 
the symmetry breaking external field is specified by 1/ cosh(iV r i?) hence it becomes smaller as the constituent quark 
mass gets larger. Thus Satz's argument on the relation between the Polyakov loop behavior and the chiral dynamics 
is already embodied qualitatively in the Gocksch-Ogilvie model in a natural way. Of course, the model will give 
us deeper insight into the interplay between the deconfinement transition and the chiral phase transition beyond 
qualitative arguments, as we will see in the subsequent sections. 

Now we shall remark about the applicable range in which the Gocksch-Ogilvie model can make sense. The limit 
of the temperature going to zero would render the model illegitimate because separate treatments for the spatial 
directions and for the thermal direction as illustrated in Fig. ^ can be allowed safely only at high temperature. 
In order to take the limit of zero temperature properly in the Gocksch-Ogilvie model, the group integration with 
respect to the Polyakov loop, as performed in Ref. 17], is indispensable, which needs somewhat complicated algebraic 
procedures. The interested temperature here is, as we will see later, around 200 MeV which corresponds to N T ~ 2 
in the present analysis. Thus we believe that it is well-founded to deal with the thermal direction separately in the 
construction of the model. The opposite limit of the temperature getting extremely high would make the quantitative 
argument meaningless simply because the strong coupling expansion does not work in a perturbative regime at such 
high temperature. The Stefan-Boltzmann law, for example, cannot be derived within the framework of the strong 
coupling expansion on the lattice. In any case, as far as the region around the critical temperature is concerned, we 
can presume that the Gocksch-Ogilvie model should be applicable. 



III. ANALYTIC CALCULATIONS 



A. Mean-Field Approximation 



The dynamics of the Polyakov loop can be understood in analogy to spin systems with the center symmetry. As 
usually performed in the calculation in the spin system, we will adopt the mean field (Weiss) approximation for the 
Polyakov loop. In this approximation, we take account of the fluctuation of the individual Polyakov loop surrounded 
by a constant mean field around it. The mean field action is taken here as 

Smf[L] = -~ ^{Tr c L(n) + Tr c tf(n)}, (5) 

n 

where x is the mean field for the Polyakov loop and treated as a variational parameter to be determined afterward 
from the extremal condition on the free energy. As for the chiral order parameter, A is simply regarded as a constant 
mean field denoted by A. The mean field free energy is then given by [26j 

= -N-«-Q({-SW[L,\] + S ml [L]) mf + ]n J VLe- s ^) 

= PUuon(x) + PfjfJ^&X). (6) 

Here (■ ■ -) m f means the expectation value with the weight specified by the mean field action. N^^ 1 " 1 is the spatial 
volume. Roughly speaking, the first term (if?5j) m f corresponds to the internal energy tending to make the system 
ordered and remaining terms are the entropy tending to make the system disordered. 
The gluonic part of the free energy can be written as 

0Uuon(x) = -2{d-l)j(J^\nI{x) S ) + x 2 -^(^\nl{x^j, (7) 

where I(x) is defined by [27j 

I(x) = JdL exp |Tr c (L + tf) = ]T det I m _ i+j (x) (8) 
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with I n {x) representing the modified Bessel function of the first kind. The quark part of the free energy in the r 
representation is written for Nf = 2 as 



At) , r, _ dim(r)iV T - 2 dim(r) 

f '"* lI,AJ " 2 2 -L— v-r-,j 2 /(a . 



r , \ i 1 jWfx;cosh(JV r J5)) 
ln[cosh(7V T ^)] - ~ — 



(9) 



I^'\x] a) is given in the fundamental and in the adjoint representation respectively by 

No oo 

n7 



c a— 1 m— — oo 

/" 27r d</> , r cos 0i 

J„(a:;a!) = / — In H cos(r 

io 2f L a i 



X COS 



and 



N, 



Me 



j JV C oo 

ei i---iN c e jl---jN c Im~i 1 +j 1 {%) ' ' ' I m ^i b+ j b (x; a) ■ ■ ■ Im-i Nc +j Nc { x ) 

a^b=l m= — oo 

+ (7V c -l)ln[l + -l/(x), 
L a J 

2i ii j , 



(10) 



„ (27T) 

^ln 
2tt 



■In 



1 



COS(0 a - 6 1 



cos(m0 a + n0 6 ) e ^ cos *»+ cos ^) 



COS ( 



^ — ^ — 'P) Im + n \ 



2x cos — 
2 



(11) 



In Eqs. (|10J) and (jl 1|> the integration in terms of <p originates from the thermal excitation of fundamental quarks. Since 
the excitation of an adjoint quark consists of a fundamental quark and anti-quark pair (see Eq. ©), two integration 
variables (f> a and 4>b appear in Eq. (fTTfl . 

Now that we have the actual expressions for the mean field free energy given by Eqs. (JZJ) and (0, all we have to 
do is find the minimum of the free energy with respect to the variational parameters x and A. Before addressing the 
numerical analyses, we shall look into the analytic expressions of the free energy further in the next subsection. 



B. Chiral Symmetry Breaking in the Confined Phase 

It would be worth considering the particular cases where we can manipulate the mean field free energy in analytic 
ways. In this and the next subsections we take the chiral limit, namely, m q = and discuss what becomes of the 
chiral symmetry in some particular cases. 

An interesting limit is that the mean field for the Polyakov loop is forced to be zero by hand, i.e., x = 0. This limit 
readily corresponds to the confined phase for adjoint quarks because the Polyakov loop serves as the order parameter 
for the deconfinement transition in the system with adjoint quarks preserving the center symmetry. To be interesting, 
even for fundamental quarks, the limit of x = should indicate the confined phase. The reason is as follows; the group 
integration with respect to the Polyakov loop would pick out only color singlet combinations in the integrand. In the 
limit of x = 0, i.e., without any external source from the mean field, the expansion of the logarithmic contribution in 
terms of the Polyakov loop gives rise to only center symmetric combinations. Actually, the condition of x = leads 
to the description in the canonical ensemble with the quark triality fixed at zero (see Ref. |2^| and references therein). 
The criterion for the deconfinement transition, therefore, regains its meaning based on the center symmetry and it 
follows that the vanishing Polyakov loop signifies the confined phase even for fundamental quarks. 

For the quark part of the free energy, taking x = makes the expressions of Eqs. (|10H and simplified as 



F und -)(0;a) = NJ (0;a) - 2(-l)^ J f A r c (0; a) 



I^(0;a)=N c (N c -l)I {0:,a)- £ / a -i(0;a). 



(12) 
(13) 



6 



Then we can make use of the following integrations |39j , 

^ A± - cos^l , \a + Va 2 -l 



l (0;a) 



4(0; a) = 



2vr 



In 



1 



271 



— In 

2tt 



a 
cos ( 



In 



a 



2q 
1 



cos 



(™/>) = (Va 2 - 1 -a) r 



(n ^ 0) 



to come by the expanded form of the free energy in terms of A as |4C 



/3/i f r d °(0,A)~ (const.) 



0f^-\Q,X)^ (const.) 



d-lfN, 



-1)N T \\\ + 0(X Z ), 



d-l 



N. 



2 1 



iV T |A| + 0(A 2 ), 



(14) 

(15) 
(16) 



where [• • • ] in Eq. i|ltj|) denotes the greatest integer function which gives the greatest integer less than or equal to the 
number. To derive the above expressions, we have substituted a = cosh(N T E) ~ 1 + (d — l)iV 2 A 2 /4. 

Here we should emphasize the appearance of a linear term of |A| which stems from the expansion of \J o? — 1. Such 
a linear term always causes instability around the chiral symmetric vacuum (A = 0). In other words, we have shown 
from Eqs. i|15|) and i|16|) that the Gocksch-Ogilvie model with the mean field approximation has the property that the 
chiral symmetry is inevitably broken spontaneously in the confined phase for both fundamental and adjoint quarks. 
This property is in accordance with the general theoretical requirement. 



C. Chiral Restoration Temperature 



We can now consider the opposite situation where the Polyakov loop is almost trivial in the deconfined phase and 
the chiral restoration occurs continuously. Then, as discussed in Refs. |2ll |22| . we can acquire the chiral restoration 
temperature by the quadratic term of the free energy expanded in terms of the chiral order parameter A. Supposing 
that the mean field for the Polyakov loop takes a substantially finite value, let us assume that the Polyakov loop can 



be approximated simply by Li 



5ij and L 



(adj.) 



• 5 a b, though this approximation breaks the unitarity of the 



matrices L and i( ad J ). The free energies for fundamental and adjoint quarks are expanded up to the quadratic order 



'(£ = «, A) (const.) + 



1 - 



(d- l)N T 
4(v + l) 



A 2 + 0(A 4 ), 



a J-( ad j ) 
< J quark 



(L = v, A) ~ (const.) + 



(TV 2 - 1)JV T 



1 - 



(d- l)N T 
A{v 2 + 1) 



A 2 + 0(A 4 ). 



From these expansions, we can immediately have the chiral restoration temperatures; 

d—1 -i „r*At i d—1 i 



/j^(fund.) 



1 



7"(adj.) _ 
X 4(^2 



1) 



(17) 
(18) 

(19) 



In contrast to the lattice observation |24| . the chiral restoration temperature for adjoint quarks is not changed radically 
from that for fundamental quarks. From the analytic expressions for the critical temperature we can see the tendency 
that the chiral restoration should be hindered as the system approaches the confined phase where v ~ (Tr c L) gets 
smaller, though the above expressions are only valid for v ~ 1. 



D. Modification to the Deconflnement Temperature 

Although the angle integrations appearing in Eqs. (|10fl and Ijllf) seem too complicated to accomplish analytically, 
the expansion in terms of 1/a, which will be acceptable when the quark mass is large, enables us to write down 
I n (x;a) and I™(x;a) immediately (sec also App. [DJ). Such an expansion in the leading order of 1/a and a further 
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FIG. 2: Schematic figures to show the modification to the free energy. The left figure is for the first order transition at T — Td 
in the pure gluonic case. The transition temperature is lowered from Td to Tj* and the minimum around the origin is shifted 
by the linear term for fundamental quarks (middle figure). The transition temperature is also lowered to T^* by the quadratic 
term for adjoint quarks (right figure). 



expansion in terms of x as well yield the following expressions of the free energy for N c — 3; 

NT /x x 2 5x 4 x 5 \ 2 x 2 x 3 x 5 5x 6 , . 

W*) - -2(d - 1) J( - + - - — - — ) 4- T + T2 - Vq - ^ (20) 



x 2 


5x 4 


8 


~ 384 


X 3 




12 


~ 96/ 



/3/ (fun ^ (x, A) ~ -N T \ 2 - - In cosh( A T T) \ (- + — - — -—), (21) 

quark V > ; 2 r 2 v r ; 2 cosh( V T T) V 2 8 384 256/ 



AfS (z, A) =s 4A r A 2 - 4 In cosh( A r £) ? f — + — - — V (22) 

MJ qua rkV . ) r \ r ) 2 COsh( V T T) V 4 12 96/ V ' 

The gluonic part of the free energy always has a minimum around x = because the expansion starts from the order 
of x 2 . Thus it can describe the first order deconfinement transition from x = to some finite x at the deconfincmcnt 
temperature Td as shown schematically in the left figure of Fig. The absence of the linear term of x reflects the 
center symmetric property. 

The expansion l|21|) for fundamental quarks clearly shows that the quark contribution breaks the center symmetry 
through the linear term of x. In the presence of such a linear term, the minimum around x = in the gluonic free 
energy is shifted to a finite value of x. As far as only the linear term, which is negative now, is concerned, the free 
energy will be more decreased at larger x. As a result, the first order deconfinement transition temperature is lowered 
from Td to Tj* < Ti by the effect of the explicit breaking of the center symmetry, if the first order transition persists. 
The intuitive picture is depicted in the middle figure of Fig. [21 As a matter of fact, the first order transition occurs 
even in the presence of dynamical quarks with m q > 300 MeV for Parameter I and m q > 700 MeV for Parameter II 
(the meaning of the model parameters is given in the next section). 

In Eq. Q22p. on the other hand, the leading term from the logarithmic contribution is quadratic in x and the 
minimum around x = is not altered then. This is because of the center symmetry remaining for adjoint quarks. 
When we focus attention on the leading term, as in the above case for fundamental quarks, the negative contribution 
of order x 2 in the quark free energy will reduce the free energy more for larger x. Therefore, as shown in the right 
figure of Fig. [21 the first order deconfincmcnt transition temperature is also lowered to T^* < Td in the presence of 
adjoint quarks. 

The strength of how much the deconfinement temperature is lowered, i.e., Ti — T^ or Ti — Tf*, is specified by 
~ 1/ coah(N T E). As long as the current quark mass, m q , is smaller than the scalar condensate A, or as long as the 
constituent quark mass Ea~ x is smaller than the temperature T = l/N T a, the modification to the deconfinement 
temperature will be almost insensitive to m q . This anticipation, in fact, has been observed in the lattice aQCD 
simulations l24l and will be confirmed in our numerical results in the next section. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



The parameters inherent in the Gocksch-Ogilvie model are the lattice spacing (cut-off) a and the current quark 
mass m q . In our previous work 22], a and m q are chosen to fit the masses of pions and p mesons. If m q is much 
smaller than a -1 then a is essentially determined by the mass difference between pions and p mesons and m q is fixed 
by the condition to reproduce the pion mass. The numerical values of the model parameters in this case are listed as 
Parameter I in TableQ] The deconfinement transition in a pure gluonic system occurs at 2(d — 1) J c = 0.806, that is; 
Ti = 208 MeV with these values of the model parameters. This value of seems to be rather low as compared with 
the empirical value (Ti = 270 MeV) found in the lattice simulations. Since the transition temperatures are crucial 
rather than the hadron spectrum here, we will try another parameter set, Parameter II, in this paper; instead of 
using the p meson mass, we can fix the lattice spacing a so as to reproduce the empirical value of the deconfinement 
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a 1 m q 


m p 




T d 


T c 


Parameter I 


432 MeV 5.7 MeV *140MeV 


*770 MeV 


-(342 MeV) 3 


208 MeV r 


- 187 MeV 


Parameter II 


333 MeV 7.4 MeV * 140 MeV 


597 MeV 


-(263 MeV) 3 


* 270 MeV r 


- 230 MeV 



TABLE I: Model parameters employed in the numerical calculations and the resultant outputs. *-quantities are used as inputs 
in order to fix a and m„. \& is the chiral condensate. 



Parameter I 



Parameter II 







<Tr c L>/A' c 


// 


pure gluonic case 




)\\ 

/' Y 
/' \ 


\ 

N 




// V 

// » 
// 













100 



200 

Temperature [MeV] 



300 



<Tr,L>/A', 



pure gluonic case 




100 



200 

Temperature [MeV] 



300 



FIG. 3: Behaviors of the order parameters as functions of the temperature for N c 
in the left figure (Parameter I) and for m q = 7.4 MeV in the right figure (Parameter II). The thin solid curve is for m q — 1 MeV 
and the dashed curve for m q — 50 MeV. 



transition temperature by substituting J c and the empirical Xd for Eq. We summarize our choice of the model 
parameters and the resultant outputs in Table [I] As we will discuss later, these two parameter sets (Parameter I and 
Parameter II) happen to correspond to two typical cases for the realization of the simultaneous phase transitions. 

It should be noted that we have determined the model parameters by using the expressions of the hadron masses at 
zero temperature according to Ref. , instead of those achieved in the zero temperature limit of the Gocksch-Ogilvie 
model. We briefly review the necessary expressions at zero temperature with modifications for Nf flavors in App. [U] 
In order to take the zero temperature limit of the Gocksch-Ogilvie model we have to integrate out the Polyakov loop 
variable properly. It is, however, a somewhat subtle procedure and thus we choose to use the would-be solutions 
at zero temperature from the beginning. Actually the numerical value of the model parameters fixed in this way 
turns out to be almost the same as that fixed by the approximate zero temperature limit (without the Polyakov loop 
integration) of the Gocksch-Ogilvie model as actually adopted in Ref. [2l| . 

The numerical results for 7V C = 3 are presented in Fig.[3|together with the results for several different choices of the 
current quark mass. In the presence of quark contributions, the expectation value of the Polyakov loop is not equal 
to x itself, though in the pure gluonic case the extremal condition on the free energy ensures that the mean field for 
the Polyakov loop is simply provided by x. From Eqs. JSJ and (JSJ we can have the following relation; 



(Tr c i- 



Tr r Lt) = 



1 dl(x) 
N~ c dx 



(23) 



from which we have calculated the expectation value of the Polyakov loop. As for the chiral order parameter, we 
simply use the scalar condensate A, for we can show that it is proportional to the conventional order parameter, i.e., 
the chiral condensate ^ = (xx)- Further details are given in Add. IbI 

In both cases of Parameter I and Parameter II, we can see the order parameters for deconfinement and chiral 
restoration indicating crossover behaviors almost simultaneously. Although the numerical results in Fig. might seem 
just supportive for Satz's argument, the relation between the Polyakov loop and the chiral order parameter has richer 
contents in itself, which can be deduced from Fig. 



9 



Parameter I 



5 
(75 



0.03 
0.02 
0.01 


0.03 
0.02 
0.01 



Ill 

Xt \ 

J r^\_ X 


m ? =5.7MeV |_ 

X 






0.03 










m,=50MeV |^ 






L 

Xt 








0.02 








^X 


X 


0.01 

















150 




200 




250 



0.02 



0.01 



Parameter II 




S 0.01 





0.02 



0.01 



1 1 

X 


m 9 =lMeV | 


A X{ 






L 




Xt 






x t \ 


m ( =50MeV | 

L 

Xt 





Temperature [MeV] 



200 250 

Temperature [MeVl 



FIG. 4: Slopes (temperature derivatives) of the order parameters. 



In Fig. 21 we plot the temperature derivative of the order parameters defined by 

1 d(Tr c L) 



Xt 



N c dT 



xt 



dT 1 



(24) 



which are nothing but the slope of the curves drawn in Fig. 

In both cases of the parameter sets, the pseudo-critical temperature, T c , inferred from the peak of Xt 1S almost 
identical to that from the peak of Xt ■ Surely the coincidence of the pseudo-critical temperatures found in the lattice 
QCD simulation can be reproduced in the slope of the order parameters within the present model study. It is 
interesting to see that the physical implication for the case of Parameter I is different from that for the case of 
Parameter II. 

For m q — 1 MeV with Parameter I, two peaks emerge apparently in the behavior of Xt ■ ^ n this case one peak around 
T c ~ 187 MeV corresponds to the remnant of the first order deconfinement transition rather than the chiral phase 
transition as conjectured formerly in Ref. |2^. The other around T ~ 210 MeV has the origin from the second order 
chiral phase transition, which can be also confirmed by the thin curve behavior in Fig. [3] Therefore the coincidence 
of the pseudo-critical temperatures signifies not the chiral phase transition but the induced instability of the chiral 
order triggered by the deconfinement transition. 

For m q — 1 MeV with Parameter II, on the other hand, there appears only one peak around T c ~ 230 MeV. The 
peak height of Xt around T c becomes larger as the current quark mass is lowered, which is peculiar to the genuine 
phase transition. Therefore, this peak is considered as resulting from the second order chiral phase transition and the 
transition temperature, in fact, agrees with that read from the thin curve behavior in Fig. [3] Thus the coincidence of 
the pseudo-critical temperatures in this case can be understood according to Satz's argument. 

Which transition governs the simultaneous jump around the pseudo-critical temperature relies on the magnitude 
relation between the chiral restoration temperature T x and the lowered deconfinement temperature T£. In the case 
of Parameter I, < T x is realized to lead to, so to speak, the deconfinement dominance around T c . In reality, 
because two peaks are not observed in the lattice QCD simulation, a more relevant situation is to be described by 
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FIG. 5: Susceptibilities with respect to the order parameters. 



Parameter II where > T x or ~ T x would hold. The former > T x means the chiral dominance around T c 
and then the critical phenomena are to be characterized by the critical exponents associated with the chiral phase 
transition. However, because of the general property that the chiral phase transition must take place at higher 
temperature than the deconfmement transition docs, as we have explained in the previous section, the latter ~ T x 
is likely to be realized. If it is the case, the chiral phase transition might be so affected by the remnant of the first 
order deconfmement transition that the critical phenomena peculiar to the second order transition would be obscured. 
This can be one possible account for the unsatisfactory agreement between the anticipation based on the universality 
argument and the critical exponents measured in the lattice QCD simulations [Tsl ]. 

Fig. El is the figure for the susceptibility with respect to the order parameters. The susceptibilities are here to be 
calculated by 



Xh 



^(<(Tr c L) 2 ) - <Tr c L) 2 ) = ^5 ^ !*/(*), (25) 

xx = (A 2 ) - (A) 2 = ( d j . (26) 

The definition of \\ seems to be different from the conventional chiral susceptibility = (ixx) 2 ) ~ (xx) 2 ''■> but \\ is 
almost proportional to x* near the critical temperature as explained in App. [bJ 

The genuine susceptibilities appear to be substantially different from the slopes and do not indicate any clear coin- 
cidence of the location of peaks. This is partly because of the very difference between the slope and the susceptibility 
and partly because of the approximation. 

The slope \t an d the susceptibility \\ are certainly sensitive to the second order chiral phase transition, and 
yet they can have different behaviors. In the chiral limit for a clear example, the slope vanishes above the critical 
temperature T x , while the susceptibility becomes finite or even larger than the magnitude below T x . Actually simple 
mean field arguments typically give twice larger amplitude for the susceptibility above T x than that below T x . 

The slope and the susceptibility relevant to the Polyakov loop, on the other hand, do not have to enhance due to 
the first order deconfmement transition. Therefore, for Parameter I, we have no reason to expect any coincidence 
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FIG. 6: Behaviors of the order parameters for adjoint quarks with Parameter I. The thick solid curve is for m q — 5.7 MeV and 
the thin solid and the dashed curves for m q — 1 MeV and for m q = 50 MeV respectively. 

from the beginning because of the deconfinement dominance where the Polyakov loop behavior has little information 
on the second order transition. The Polyakov loop behavior is still smooth even for Parameter II. We consider that 
it is because the simple ansatz of the mean field action {JJ would miss the significant part of the coupling to the 
symmetry breaking contribution through which the coincidence is brought about. This ansatz can incorporate, at 
most, a linear term of an external magnetic-like field and thus the non-linearity leading to the fact that the chiral 
phase transition temperature must be higher than the deconfinement transition temperature would be lost as a result 
of the approximation. Hence we cannot see the coincidence of the peak location even for Parameter II. 

Now we would comment upon one interesting possibility: The results from the lattice QCD simulation suggest that 
the situation of Parameter II is more realistic, and nevertheless, it would be appealing to consider some modification 
or artificial tuning in the lattice simulation, by which the case of Parameter I is realized. Then we can make use of 
the lattice simulation directly to investigate the nature of the deconfinement transition around T c and the nature of 
the chiral phase transition around T x > T c separately. 

Finally we show the numerical results for adjoint quarks with Parameter I in Fig. E] The calculation has been done 
with the function I™(x; a) expanded in terms oil/ a up to third order, for the numerical integration is not impossible 
but somewhat tough. We give the explicit expressions and verify the reliability of such an expansion in App. [D] As 
expected, the deconfinement transition of first order persists clearly and the chiral phase transition occurs distinctly 
at higher temperature than the deconfinement transition does. 

The global features consist well with the lattice aQCD results. The deconfinement temperature is lowered from 
Td = 208 MeV to ~ 180 MeV. For various choices of m q , the lowered deconfinement temperature T^* is almost 
the identical, as we have already discussed in the previous section. 

The behavior of the chiral order parameter, as pointed out in Ref. [24j], changes radically at T — the current 
quark mass dependence is not manifested until the temperature exceeds the deconfinement temperature. This feature 
can be understood qualitatively from Eq. I|19|) . While the Polyakov loop is zero, the chiral order parameter behaves as 
a function of the temperature as if the chiral phase transition would take place at much higher temperature. After the 
deconfinement transition, the chiral restoration comes to affect the chiral order parameter behavior and, as a result, 
the quark mass dependence becomes apparent then. 

V. CONCLUSIONS 

We investigated the relation between the deconfinement transition and the chiral phase transition at finite tem- 
perature in the Gocksch-Ogilvie model with fundamental or adjoint quarks incorporated. We shall enumerate our 
findings here; 

1. We have shown in an analytic way that, in the confined phase with the vanishing Polyakov loop, the chiral sym- 
metry should be spontaneously broken at any temperature in the Gocksch-Ogilvie model for both fundamental 
and adjoint quarks. 
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2. We discussed qualitatively how the quark contribution to the effective action can lower the deconfinement 
transition temperature. 

3. In the case of fundamental quarks we have acquired the behavior of the order parameters as functions of the 
temperature and their slopes by differentiating with respect to the temperature. Then we found the peaks of the 
slope for the Polyakov loop and for the chiral order parameter located at the same pseudo-critical temperature. 

4. Physical implication of the coincidence of the pseudo-critical temperatures has proved to be twofold: One is 
the deconfmement dominance where the remnant of the first oder deconfinement transition would govern the 
behavior of the order parameters around T c . The other is the chiral dominance where the second order chiral 
phase transition would be responsible for the realization of the simultaneous jump of the order parameters around 
T c . We speculate that the real world corresponds to an intermediate admixture between the deconfinement 
dominance and the chiral dominance. 

5. In the case of adjoint quarks we have found that the deconfinement transition is seen clearly of first order and 
that the chiral phase transition of second order occurs at higher temperature than the deconfinement transition 
does. It is apparent that the deconfinement temperature is substantially lowered from to T^*, as expected, 
and the difference Td — T^* is insensitive to the current quark mass m q as far as m q is smaller than the scalar 
condensate (i.e. the constituent quark mass) or the temperature. 

Some comments arc in order. We will address our conjecture here that the finding 1 would be generalized so that 
we can prove the chiral symmetry breaking in the confined phase only using the property of the group integration, 
or, the Haar measure. Also it would be interesting to make sure whether 1 is realized in the lattice simulation. Such 
simulations will be feasible, in principle, if once we can examine the chiral symmetry breaking with a constraint 
imposed on the expectation value of the Polyakov loop, i.e., (Tr c L) = 0. 

As for 3, we must improve the mean field approximation, especially the ansatz of the mean field action, to reproduce 
the coincidence of the pseudo-critical temperatures in the genuine susceptibility. This is a task for the future. 

It would be fascinating from 4 to consider the possibility of the deconfinement dominance realizing in the lattice 
simulation because we can investigate the nature of the deconfinement around T c ~ T£ and the nature of the chiral 
restoration around T x > T c separately. In order to achieve such a situation, we have to make some modifications to 
the lattice QCD either to lower or to raise T x . The problem of how to put the idea into practice is beyond the 
scope of the present paper. 

Although it is rather academic, we finally remark about the possible relevance to physics on the critical end point 
(CEP). Recently the CEP, that is the terminal point of the first order transition, is paid much attention to in the 
realistic QCD phase diagram l28l. The location of the CEP is determined in part in the lattice QCD simulation at 
finite temperature and density |29j . The CEP can also emerge as the strange quark mass is varied in the three-flavor 
chiral phase transition |3Cj . 

Here we will draw attention to another chance to be confronted with the CEP, namely, the point where the 
deconfinement transition ceases to be first order with massive dynamical quarks. Our numerical analysis in the 
Gocksch-Ogilvie model predicts that the critical quark mass m q ~ 300 MeV for Parameter I and m q ~ 700 MeV for 
Parameter II. Actually there exist some lattice measurements to determine the critical mass of the CEP: In Ref. [3l| 
m q ~ 900 MeV is concluded, while Ref. j32 reads m q ~ 400 MeV in more improved analyses with two-flavor quarks. 
The results in Ref. [12 is consistent with our estimate of m c q . In the one- flavor QCD simulation m q ~ 1.4 GeV 
is reported |3^|. As for the three-flavor case [34[, the situation becomes somewhat subtle because not only the 
deconfinement transition with heavy quarks but also the chiral phase transition with light quarks can be of first order. 
We anticipate that the property of the CEP in the deconfinement transition might be intriguing from the point of 
view of the glueball physics at finite temperature [3^, |3(| . 

In view of the prosperous descriptions as listed above, we can regard the Gocksch-Ogilvie model as a promising 
implement to afford a conceptual clue to the relation between the deconfinement transition and the chiral phase 
transition. It would be intriguing to make progress in applying the model to the system at finite density. 
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APPENDIX A: CONSTRUCTION OF THE GOCKSCH-OGILVIE MODEL 



In order to make this paper self-contained, we shall explain how to construct the Gocksch-Ogilvie model explicitly 
by following the derivation of Refs. fl6l . l2l| . This would be useful to clarify what approximation is implicitly presumed. 
The model is formulated on the lattice. The strong coupling expansion is readily feasible on the lattice to lead to 
the quark confining property, i.e., the area law of the Wilson loop. Although the real world is actually in the weak 
coupling regime, numerical studies in the lattice QCD have revealed that even the strong coupling expansion should 
be likely to give a plausible picture relevant to the real world. The spontaneous breaking of the chiral symmetry 
can be described on the lattice as well by means of the double expansion of strong coupling and large dimensionality 
[ilj Il7j . Thus we can expect that an effective action suitable for discussing the relation between the Polyakov loop 
and the chiral order parameter will be available if once we bring together those assets on the lattice. This is the 
starting point to construct the Gocksch-Ogilvie model. 

We adopt the staggered fermion in order to approach the chiral limit in a plain manner and also in order to take 
advantage of the simplicity in regard to the spin structure. The action on the lattice is given by the gluonic part; 



S G [U] 



n,(/i.,i^) 



1 



i 



ReTr c £/ M1/ (n) L U^(n) = Ufa)Ufa + i>)U v (n + fiU^n) 



(Al) 



for N c > 3 and the staggered quark part; 



s F [u d ,Ui,x,x] = m ?y^x(«-)xW + 



- E 

2 ^ 

n,/x>0 

l x 



, n,x{ri)x{n){U il (n)x{n + A) - Ufa - p) X {n - p.)} 



= m q E X(n)x(n) + 7; E Vd(n)x(n){U d (n)x(n + d) - Ufa - d)x(n - d)} 

n n 

Tr c { $i (n)Ui (n) + $i (n) U} (n) } , 



(A2) 



n,i^d 



where 



{^( n ))ah = -\vi(n)Xb(n)xa(n + i). 



(A3) 



We have divided the quark action into one part with Ui and the other part without Ui for later convenience, a, b run 
in the color space. We can handle the quark action in the adjoint representation by replacing the link variables in the 
covariant derivatives by 



U, 



(adj.) 



ab 



2tr 



t a Ut b tf 



(A4) 



Roughly speaking, the adjoint quark is like a color octet fermionic meson composed of a quark and an anti-quark in 
the fundamental representation (8 out of 303 = 108). We will write the link variable as in the r representation. 

In order to attain the effective action in terms of the Polyakov loop and the quark field, we must integrate over the 
spatial link variable Ui. Then the effective partition function is written as 



z^[u d , x ,x} = 



exp 



X(n)x(n 



Vd{n)x(n){U { d r \n)x{n + d) - U d r) \n - d) X (n - d)} 



VUi exp 



S G [U d ,Ui] + TT c [Mn)U^(n) + ^(n)U^\n)} 



(A5) 



In the strong coupling expansion, the spatial plaquette can be neglected, which can be perceived more on the 
asymmetric lattice |37| . An extreme case of such asymmetric lattice methods is the Hamiltonian formalism where 
only the temporal plaquette remains in the strong coupling limit. Then the Ui integration part can be expanded as 



T 1 p ' 

f VU * E i 3 E ^{U id (n) + Uj d (n)\ £ Tr c {^ r \n) + 9^)17™ (n)\ 



(A6) 
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FIG. 7: Graphical representation of the configuration contributing to the effective action in the lowest order of the strong 
coupling. The torus stands for the thermal (temporal) S 1 and the spatial R 3 . 



The Taylor expansion of the gluonic part (series in p) obviously corresponds to the strong coupling (l/g 2 ) expansion. 
As we explain soon below, the Taylor expansion of the quark part (series in q) corresponds to the large dimensional 
{1/d) expansion. We specify the order of the double expansion as follows; first, we perform the large dimensional 
expansion. Then, in each expanded term we take the leading contribution of the strong coupling expansion. 

For the zeroth order (q = 0) term of the 1/d expansion, the leading contribution left after the Ui integration derives 
from the configuration with the plaquettes wrapping around the thermal (temporal) torus (p = N T ), as depicted in 
Fig. |7| As a result, the Ui integration for the q = term is ^3 



1 



v5 2 ^c 



^Tr c L(n)Tr c L\m). 



(A7) 



The next contribution comes from the second order (q = 2) term. Since two link variables in the covariant derivative 
can form a color singlet pair, the leading contribution in the strong coupling expansion is simply the zeroth (p = 0) 
term. As a result, we have 



JvU z Tr {*iU$ r \n) + $,(n)t/f r)t (n)} 

\_n,i^d 



2(d-l) 



M(n)M(n + i), 



(A8) 



where the mesonic composite field is defined by 



M(n) = 



dim(r) 



— dim(r) 



a=l 



Xa{n)xa{n). 



(A9) 



In Eq. I|A8(1 the dimension number of the r representation is denoted by dim(r) as in the text. It should be noted that, 
for adjoint quarks, only such a mesonic interaction as Tr c {$j$j} survives not only because of the group integration 
with respect to Ui but also because of the properties of the Grassmann variables x an d X- 

At this stage it is easy to understand why such an expansion in q provides the large dimensional expansion. Since 
the meson should propagate with finite amplitude in the large dimensional limit, we must properly normalize M(n) 
so as to have the factor l/(ef— 1) in front of the term M(n)M(n + i) (i runs from 1 to d — 1). As a result, each 
time $i, or M, appears in a term, the order of the term in the 1/d expansion should be lowered. Thus the Taylor 
expansion of the exponential of the quark action amounts to the expansion in 1 / d. 

We combine Eqs. I|A5|1 — (|A8|) and exponentiate the expanded terms. Then, using the identity, 



exp 



dim(r 



T>\ exp 



-dim(r) | — A(n)V(n, m)X(m) 



X(n)V(n, m)M(m) 



(A10) 
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we can reach the effective partition function; 



Z^[U dl \] = J V X Vxexp 



S (G) [L] dim(r) ^ x{n)v ^ m)x{m) 

m,n 

- ~ X>(n){x(n)E#V)x(« + d ") " X(" + <2)Etf° + (n)x(n)} 

n 

- Z j y Z Hm)V{m, n) + m\ X {n d ) . (All) 

n \ m ) 

From Eq. (| AT|> . the effective action in terms of the Polyakov loop, S^[L], is given by 

S$[L] = -J^Tr CJ L(n)Tr c Lt(^). (A12) 

n.n. 

where J is given by 

J = f4rf T = e-' o/r . (A13) 



Here we have used the familiar expression for the string tension in the strong coupling expansion, i.e., era 2 = ln.[g 2 N c ]. 
We will fix the lattice spacing a as a cut-off parameter and vary the temperature by changing the number of lattice 
sites in the thermal direction, N T . 

In Eq. (|A11|) the hopping propagator of the meson fields V(n,m) is defined in Eq. Since the quark part takes 
a bilinear form in Eq. JIllJ), we can perform the integration with respect to the quark fields explicitly under the 
approximation that A is independent of rid, (time) . This approximation is indeed reasonable as long as we employ the 
mean field approximation and treat A as a constant, in the same sense as the sigma meson field in the relativistic 
nuclear many-body system. We would emphasize, however, that this approximation might be invalid when the meson 
fluctuation is concerned. In other words, we should not use the action of the Gocksch-Ogilvie model to estimate the 
meson fluctuation (loop) effect. 

Then after the integration we can reach the Gocksch-Ogilvie model given in Eq. (JJJ by using S^g = — In . 
We should remark that N{ in the effective action Q is the number of flavors introduced from the phenomenological 
point of view. We cannot avoid having four degenerate flavors per a single component of the staggered fermion. The 
logarithmic contribution in Eq. (JJJ results from the quark integration in Eq. (|A1 lfl . The coefficient iVf/4 in front of 
the logarithmic term effectively adjusts the number of flavors in the integration with respect to the quark fields x- 

It is worth mentioning that a in Eq. (|A13|) is constant and fixed at the empirical value at zero temperature. The 
string tension of the linear potential defined by the Polyakov loop correlation function, nevertheless, depends on the 
temperature and decreases with the temperature increasing. In reality the string tension should have temperature 
dependence in addition due to the vibrations of the flux tube. The higher temperature the system is heated up at, 
the more furious the vibration becomes. Such effects of the string vibration are to be described by the higher order 
contributions in the strong coupling expansion. If we take account of the spatial link variables, or the magnetic 
plaquettes in the Hamiltonian formalism, to exceed the leading order contribution of Eq. I|A7|) . we can incorporate 
temperature dependence coming from the string vibration in principle. Since the essence on the deconfinement 
transition is already concentrated in the effective action (|A7|) . we simply neglect any spatial fluctuation in the present 
model. 

If we proceeded to the higher order contribution of the strong coupling expansion, we would have the fat meson 
contributions as shown in Fig. |S| Neglecting such fat excitations implies that we neglect the structural change of 
mesons in the temporal direction. In other words, the meson hopping in the spatial direction have no extent in the 
temporal direction and is always color-confined. 



APPENDIX B: CHIRAL CONDENSATE 



The chiral condensate can be calculated as 

dim(r)iV d dm q cff dim(r) dm q ' ^ ^ 
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FIG. 8: Examples of neglected excitations as a result of the approximation. The vertical and horizontal lines represent the 
temporal (thermal) and spatial link variables respectively. 



where Ao is the solution of the stationary condition (see Eq. ©) 

= dim(r) - 



A) 



dX 



\=X 



1 df^(x,X) 



di 



(B2) 



A=A 



Then, using the stationary condition, we can readily obtain 

1 d/«(x,A ) _ 1 df^{xX) 



dim(r) 
1 

dim(r) 



dm n 



dim(7') 
d\ 



dm„ 



dim(r)Ao 



-A . 



(B3) 



Actually this result does not depend on the concrete form of the free energy. We can immediately reach the same result 
from Eq. I|A10|) by infmitesimally shifting the integration variable. Thus, the scalar condensate is simply proportional 
to the chiral condensate. 

In the same way, we can relate our susceptibility \\ defined in Eq. (|26(l to the conventional chiral susceptibility 
which is given by 



1 



X* = 



[dim(r)] 2 
1 



dm 2 



[dim(r)] 2 I i dm q 



d\ 2 



mf ' 

dm 2 



(B4) 



Differentiating the first equation of Eq. (|B2J| by A and m q respectively and combining them, we can reach the relation; 



dA n 



dim(r} 



d-1 



OX 2 



d-1 



(B5) 



In the vicinity of the critical temperature with small current quark mass (m q ~ 0), the curvature of the free energy 
(second partial derivative of the free energy) with respect to the scalar meson field approaches zero. Consequently 
the first term in the brackets in Eq. 1B4JI dominates to result in 



X* 



d-1 



dX 2 



J2f3_ 
d- 



T ((A 2 >-(A) 2 ). 



(B6) 



APPENDIX C: ZERO TEMPERATURE EXPRESSIONS FOR N t FLAVORS 



Here we briefly summarize the zero temperature expressions given in Ref. [Tfjj with modifications for Nf flavors. 
The effective potential in terms of the chiral order parameter A is modified at zero temperature into 



VeS (X) = f X 2 
from which the stationary value of A is obtained as 

- 771 

A = 



N c N f 



In 



X + \ j -m q 



'2 
d' 



2d 



N { m 2 
4 2d' 



(CI) 



(C2) 
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If we rescale m q — m q \J Nf/4, then the Nf dependence in A can be factorized as 

(C3) 




where the expression in the brackets is nothing but that of Ref. |16| without any care about the flavor number. Thus 
we can expect that most of modifications could be absorbed in the rescale of the current quark mass m q — * rh q . In 
fact, the hadron spectrum for A/f flavors is given by 



coshM, = di^(^L+\l— + -^ | !K 2 P + I 



7?+ 


ra 2 q 




~2d 




2 








) 







= ^ -7= + \/l + T7j) -1K2P+1- (04) 



Therefore the hadron spectrum necessary for fixing the model parameters is not affected by the number of flavors. 

APPENDIX D: EXPANSIONS 

We shall discuss the expansion in terms of 1/a for the integrations I n (x;a) and I™(x;a) given by Eqs. (|10f> and 
(|ll|l . Such expansions are exploited for the numerical calculation with adjoint quarks. It is somewhat inefficient to 
attack the integration given in Eq. I|ll|) directly and thus we expand the free energy in terms of 1/a. We can obtain 
such an expansion from Eqs. (|10|) and l|ll|) up to third order as 

I n (x;a) = -^—{l n+1 {x) + J„_i(x)} - -^{l n+2 (x) + 2I n (x) + I n - 2 (x)} 
2a 8cr 

+ TTT^I^+sO) + 3I n +i{x) + 3I n -i(x) + I n _ 3 (x)}, (Dl) 

I™(x;a) = — {l m+ i(x)I n -i(x) + I m -i{x)I n+1 (x)} 

- -^-^{I m+ 2{x)l n -i{x) + 2I m (x)I n (x) + I m -2{x)I n +2{x)} 

+ -^^{l m +z{x)I n - 3 {x) + M m+1 (x)I n -x(x) 

+ 3I m - 1 (x)I n+1 (x) +I m . 3 {x)I n+3 (x)}. (D2) 

We can test the reliability of such approximations by comparing with the full (not expanded) numerical results. In 
Fig. we show the full numerical result obtained from the free energy © directly (dashed curve) and the approximate 
one with the expanded free energy (solid curve) in the case of fundamental quarks with Parameter I. It is obvious 
from the figure that the expansion I|D1|) works well except for the behavior of the chiral order parameter around the 
chiral phase transition. Thus we can presume that the expansion l|D2|l will work enough for adjoint quarks as well, 
at least, to see the qualitative features. 
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